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Abstract 
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energy surface of Aug. These structures provide new insights on the 2D =^ 3D tran- 
sition in small neutral and charged gold clusters. It is demonstrated that the size 
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the 2D - 3D coexistence develops for Au^ and Au^ on the cationic potential energy 
surfaces while only for Aug on the neutral. Two metastable long-lived dianions of gold 
clusters are also reported. 



^Address for correspondence: FAX: +32 (4) 366 3413; E-mail addresses: fremacle@ulg.ac.be and 
eugene.kryachko@ulg.ac.be. 
^Maitre de Recherche, FNRS (Belgium). 



1 



I. Introduction 

Finite gold clusters are very common building blocks for nanostructured materials, 
electronic devices and novel nanocatalytic systems and as such, they have triggered 
numerous experimental and theoretical studies (Refs. [1, 2] and references therein). 
Some of them have shown that neutral and charged gold clusters favor 2D structures 
to unusually large sizes [3, 4] due very strong relativistic effects [5]. The motivation of 
the present work is to investigate in details the 2D - 3D coexistence for neutral clusters 
of gold and to extend it on the positively and negatively charged ones. Before discussing 
our own results and the new structures that we identify, we start this Introduction with 
a brief review of previous works on small gold clusters. 

In one of the first ab initio calculations of small neutral and charged clusters of 
gold Au4,5, based on the size-consistent, SCF-based modified coupled pair functional 
(MCPF) method in conjunction with the Hay-Wadt 19-(5s^5p^5(i^°6s) valence electron 
relativistic effective core potential (RECP), Bauschhcher et al. [6] predicted the most 
stable Au5 cluster to be a planar W-shape of C2v symmetry (see Figure 1). This 
geometry was later confirmed by Bravo-Perez et al. [7] within the second- and fourth- 
order perturbation M0ller-Plesset (MP2 and MP4) method with the same Hay-Wadt 
RECP and by Gronbeck and Andreoni [8] who used a variety of density functional 
(DF) methods. The 3D Aus clusters were reported to be higher in energy. Two low- 
energy trigonal bipyramid structures were found 0.42 and 0.95 eV above the W-shape 
in Ref. [7] (0.79 eV in Ref. [9] within the LDA PW91 method). The LSDA method 
with the Perdew-Zunger parametrization of the Ceperley-Alder data for the electron 
gas [8] also favors the trigonal pyramid as the most stable 3D structure of Aus, though 
0.72 eV higher than the W-form. The square pyramid is found 0.87 eV above the 
W-structure [8] . However, the BLYP density functional with the same RECP shghtly 
reverses this order: the square pyramid lies 1.00 eV above the W-form while the trigonal 
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pyramid is somewhat higher at 1.03 eV [8]. Using the BP86 DF method together with 
the Stuttgart 19e-RECP, Bonacic-Koutecky et al. [10] found that the trigonal 
bipyramid is 0.930 eV higher than the W-structure. For anionic five-gold clusters, the 
square pyramid becomes more stable compared to the trigonal one (0.88 vs. 1.14 eV 
for LSDA and 1.05 vs. 1.25 eV for BLYP) [8]. 

The global minimum of the six-gold cluster possesses a planar, Ds^-symmetry tri- 
angular structure [7]. The pentagonal pyramid is energetically higher by 0.44 eV [7]. 
Due to this smaller 2D - 3D energy difference compared to that predicted for five-gold 
clusters, it was then suggested that "probably at n = 7 a 3D form could be the most 
stable or at least in degeneracy with a planar stable isomer" (Ref. [7b], p. 663). How- 
ever, as reported in Ref. [10], the pentagonal pyramid hes 0.914 eV higher than the 
triangular cluster. In contrast, using the generalized gradient approximation (GGA), 
the pentagonal bipyramid is found to be the lowest energy structure, 0.16 eV below 
the planar hexagonal structure [9]. 

Hakkinen and Landman [3, 4] theoretically predicted that the 2D =^ 3D transition 
occurs for Au„ and Au~ when n > 8 and n > 7, respectively (see also Ref. [11] for 
a study of the pyramidal Auy cluster with D^h symmetry and Ref. [12] for recent 
studies of the adsorption of methanol on gold clusters) . It has recently been shown by 
Bonacic-Koutecky et al. [10] that 3D structures appear to be higher the 2D ones for 
n = 7 - 10 by 0.528 eV (n = 7), 0.530 eV {n = 8; for opposite results for n = 7 and 
8 see the recent work [9]), 0.224 eV {n = 9), and 0.285 eV (n = 10). Interestingly, 
the differences in the binding energies per atom between the 2D ground-state and the 
lowest-energy 3D structures calculated in this work [10] fall from 0.19 eV for n = 5 
and 0.15 eV for n = 6 to 0.02 - 0.03 eV for n = 7 — 9. A simple structural rule 
has been proposed in this study [10]: the ground-state structures of Au4<„<7 can be 
derived by adding a single Au atom to the lowest-energy structure Au„_i and forming 
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an additional triangular subunit. This rule has been suggested to be a manifestation 
of the strong s — d hybridization and of the preference for the directional d-character 
bonds over the isotropic s-ones. However, this structural motif disappears for Aug 
where the square is capped by four atoms of gold on each side [10]. The 2D - 3D 
boundary for the neutrals was recently further extended to at least n — 13 [3]. The 
recent gas-phase mobility experiments by Furche et al. [13] have demonstrated that 
planar anionic gold clusters are more stable for sizes as large as n = 9 — 11. On the 
other hand, similar experiments on the cationic gold clusters, supported by the GGA 
DF computations with the S-VWN+Becke-Perdew (BP86) parametrization together 
with the Stuttgart RECP, showed that Au+ are 2D for n = 3 — 7 and 3D for n = 8 — 13 
[14]. 

Understanding the 2D - 3D transition in small neutral and charged gold clusters 
is the main theme of the present work. The paper is organized as follows. Section 
II is devoted to a brief outline of the computational approach that is based on three 
commonly used relativistic effective core potentials. In Sections III - VI we provide 
a systematic analysis of the potential energy surfaces (PESs) of Auf<„<8(Z = 0, ±1) 
and identify therein a variety of novel structures that bring new insights on the 2D =^ 
3D transition. In Section VII we demonstrate the 2D - 3D coexistence on the PES of 
neutral Aug. To conclude, the emerging trends of the 2D-3D transition as a function of 
the size, charge and odd or even number of electrons of the small gold clusters studied 
here are discussed in Section VIII. 

II. Computational Methodology 

All computations of gold clusters reported in the present work were carried out with 
the hybrid density functional B3LYP potential in conjunction with a variety of energy- 
consistent 19-(5s^5p^5(i^°6s) valence electron RECPs including the Los Alamos double- 
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zeta type RECP LANL2DZ (= A) [15], the RECP developed by Ermler, Christiansen 
and co-workers (= B) with the primitive basis set (5s5p4d) [16], and the RECP of the 
Stuttgart group (= C) with the basis set (8s7p6d) /[6s5p3d] [17] (see also Refs. [6, 7, 
18] for the description of these RECPs). The reliability of the chosen RECPs has been 
assessed in Refs. [7, 10, 13, 14, 18-22] by comparing the computed results for small 
gold clusters (dimer, trimer and their charged species) with the experimental data. 
The GAUSSIAN 03 package of quantum chemical programs [23] was used throughout. 
All geometrical optimizations were conducted with the keyword "tight" . 

The B3LYP DP has recently been used for the gold clusters [20, 24-26] along with 
the others DFs: BP86 [10, 14], BLYP [8], PW91 [19], BPW91 [27], BPW [18]. Com- 
pared to the MP2 level, B3LYP was found [20] to give very similar geometries for 
clusters containing gold atoms, i.e., the Au-Au distances are ^ A % longer than the 
MP2 ones, with the B3LYP ones being more accurate, so that the density functional 
approach might be more reahstic for the gold clusters than MP2 [20]. On the other 
hand, Bonacic-Koutecky et al. [10] and Lee et al. [18] have recently pointed out that 
the DP methods determine relatively small s — d energy gaps and underestimate the gap 
EffL between the highly occupied molecular orbital (HOMO) and the lowest unoccu- 
pied molecular orbital (LUMO). Por example, the s — d orbital energy gap between the 
and d^s"^ configurations of Au~ is estimated in the present work as equal to 0.88 
eV (B3LYP/A), 1.12 eV (B3LYP/B), 1.37 eV (B3LYP/C), 0.23 eV (BP86/A), 0.43 
eV (BP86/B), and 0.62 eV (BP86/C). These values are smaller than the experimental 
magnitude of 1.34 eV [28] except for the third estimate demonstrating a good per- 
formance of the B3LYP/C computational level for this task. To further pursue such 
comparisons for the geometrical, energetic, and MO patterns of larger gold clusters 
and to strengthen the credibility of the present conclusions, the BP86/A and PW91/A 
methods are also used for some key structures studied in this work. It is the well- 
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known fact that DF methods favor planar structures for gold clusters (see Refs. [10, 
18] and references therein). Since the present work focuses on the 2D - 3D transitions 
in the neutral and charged gold clusters studied by means of these methods, we have 
to estimate the DFT/RECP error in energy within which given Au clusters are treated 
as almost energetically equivalent. 

A value of the DFT/RECP error in energy based on the BP86/C method was 
proposed in Ref. [10]. It is defined as the energy difference between 2D and 3D 
structures for Aug clusters and equal to ^ 0.2 eV («i 4 kcal-mol~^; on p. 3124 of 
Ref. [10] this error is equal to 0.224 eV or 5.16 kcal/mol). As reported above, such 
a magnitude is in agreement with the error found in estimating of the s — d orbital 
energy gap at the B3LYP/B and B3LYP/C compared to the experimental data. We 
suggest to take this value as the error margin for the B3LYP/RECP methods used here 
to study the 2D - 3D transitions in neutral and charged gold clusters. We futhermore 
demonstrate in Section VII, devoted to the neutral Aug clusters, that the two criteria, 
viz., the energy difference between 2D and 3D structures as used in Ref. [10] and the 
s — d orbital gaps lead indeed to the same value of 0.19 - 0.2 eV (fti 4.33 - 4.55 kcal/mol, 
see Table 2). 

The harmonic vibrational frequencies and corresponding zero-point vibrational en- 
ergies (ZPVE) were calculated at all employed B3LYP/X levels (X = A, B, and C) in 
order to determine the topology of the stationary points. The electronic energies calcu- 
lated by means of the B3LYP DF with the basis sets A, B, and C are hereby reported 
in the following order: EnergyBSLYP/A (Energy bslyp/ A + ZPVE) [Energy bslyp/b 
(Ener syB3LYP/B + ZPVE); Energy bslyp/c (Energy bslyp/c + ZPVE)]. The BP86/A 
values are given in curly brackets whereas the PW91/A ones appear in the double curly 
brackets. Similar notations are used for the optimized geometrical parameters. 

The vertical electron detachment energy, VDE, is calculated as the total energy 
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difference between tlie anionic and neutral clusters, both taken in the anionic optimized 
geometry without the ZPVE. The adiabatic electron affinity, EA^ and ionization energy, 
IE(j, are calculated as the total energy difference between the anionic and cationic 
clusters, respectively, and the parent neutral cluster in their optimized geometries. 

III. Auf 

The PES of the neutral gold cluster Aus includes four low-energy planar structures 
that are shown in Figure 1. In full agreement with all the earher studies [4, 6-10, 19, 
29], at all the employed computational levels, the global minimum is the nonpolar, 
trapezoidal W-shape cluster Aug of symmetry in the electronic state ^Ai. Its 
geometry is compatible with the triangular motif rule. Au| has one 4-coordinated 
Au atom and two 3- and two 2-coordinated ones. The Au-Au distances in Aug are 
partitioned into four classes (see Table 1): the largest ones, 2.83 - 2.87, A separate 4- 
and 3-coordinated atoms; those between 3-coordinated atoms are 2.74 - 2.78 A, and 
those between 2- and 3-coordinated atoms 2.72 - 2.75 A. The 2- and 4-coordinated 
atoms are connected by the shortest bonds with the bond lengths of 2.68 - 2.70 A. The 
spin-up and spin-down gaps, ^hl-\ and ^hli, are equal to 1.58 (1.57) and 1.01 (1.02) 
eV, respectively at the BP86/A (PW91/A) computational levels (cf. with ^hl^ = 1-51 
eV in Ref. [9] and 1.37 eV in Ref. [29]). The B3LYP DP method estimates them at 
about 2.5 and 2.2 eV, respectively (Table 2). 

The next structure, by increasing energy, is the X-shape nonpolar cluster Aug^ in 
the ^B2u state of symmetry, lying 7.3 - 8.2 kcal-mol~^ above the ground W-shape 
state (see Table 2). These energies fall within the earher reported range: 7.84 - 14.53 
kcal-mol"^ [8] and 9.75 kcal-mol~^ [10]. Compared to Aug, Aug^ possesses a narrower 
spin-up HOMO-LUMO gap (0.77 eV for the BP86/A; 1.4 - 1.5 eV for the B3LYP) and 
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a slightly wider one for the spin-down electrons (1.31 eV for the BP86/A; 2.7 - 3.0 eV 
for the B3LYP; see Table 2). 

Two higher-energy clusters, Au^^^ and Au^^ (Figure 1), are reported in the present 
work for the first time. As follows from Table 2, they fill the energetic gap between 
the 2D clusters Au| and Aug^ on one hand and, on the other, the 3D pyramid which 
is ca. 21.4 kcal-mol"-*^ above Aug [10]. Both novel structures are less compact, less 
coordinated, polar (1.2 D and ~ 4.2 D, respectively) and characterized by rather large 
polarizabilities (larger by a factor of 1.5 - 1.6 compared to that of Auf). Their HOMO- 
LUMO gaps are given in Table 2. One of them, Au|^^, is rather important since it 
gives rise to the most stable anion Aul^^~ (Figure 1). Aug^^" is a novel structure too 
(see Tables 1 and 2) and is 5.8 - 7.2 kcal-mol~^ lower in energy than the well-known 
cluster Aug" [4, 8, 13, 30], which is well beyond the estimated error of ca. 4 kcal-mol"-*^ 
of the DF/RECP methods. 

The lower energy of Aug''^''^" is readily explained by its larger interatomic distances 
r(Au2-Au5) = 6.807 [6.803; 6.889] A and r(Au2-Au3) = r(Au3-Au5) = 5.300 [5.250; 
5.296] A, which reduces the mutual Coulomb repulsion between the partial excess 
electron charges that reside mainly on these atoms (the MuUiken charges are equal to 
-0.37 [-0.49; -0.55], -0.23 [-0.28; -0.23], and -0.37 [-0.49; -0.55] on Aus, Aug, and Aug, 
respectively). On the other hand, in Aug", the excess electron spreads over Au2 (-0.24 
[-0.19; -0.12]), Au3 (-0.31 [-0.40; -0.26]), Au4 (-0.31 [-0.40; -0.26]), and Aug (-0.24 [-0.19; 
-0.12]), where r(Au2-Au3), r(Au2-Au5), r(Au4-Au5) are only within 2.72 - 2.78 A. In 
Aug^^", the distance r(Aui-Au4) = 2.7 - 2.9 A is smaller by factor of two than the van 
der Waals radius of Au implying that Aui forms the Au-Au bond with Au4 and hence, 
Au5^^~ is distinct from the V-structure shown in Figure 4 of Ref. [18] (see also note 
[31]). 

These two low-energy anions Aug" and Aug^^" also differ by their MO patterns and 
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therefore have different catalytic reactivities. Within the B3LYP method, the AU5 
HOMO-1 - HOMO and HOMO-LUMO gaps are equal to 0.94 [1.01; 1.00] eV (see 
note [32] and Table 2) and 1.21 [0.93; 1.09] eV, respectively. In Au^^"" the B3LYP 
HOMO-1 - HOMO gap is substantially narrower, viz., 0.27 [0.23; 0.30] eV while the 
HOMO-LUMO one is wider and equal to 2.88 [2.39; 2.41] eV. The EA„s and VDEs of 
these anionic clusters are the following: EA„(Au5") = 3.18 [3.02; 3.07] eV; EA„(Auf 
= 4.09 [3.96; 3.99] eV; VDE(Au|-) - 3.29 [3.09; 3.16] eV and VDE(Auf = 4.20 
[4.09; 4.06] eV. The EAq and VDE of Auf~ are in agreement with the experimental 
and earlier theoretical data: EA^^p*(Au5) = 2.992 ± 0.050 eV [33], 2.93 ± 0.10 eV 
[34], 3.06(3) eV [30]; EA*„'^^'''-(Au5) = 3.33 eV [6], 3.33 eV (LSDA [8]), 3.04 eV (BLYP 
[8]), 3.06 eV [27]; VDE^^f*(Au5) = 2.98 ± 0.2 eV [33], 3.12 ± 0.05 eV [34], 3.09(3) 
eV [30]; VDE*'^^'"-(Au^) = 3.35 eV (LSDA [8]), 3.08 eV (BLYP [8]), 3.2 eV (obtained 
with the BP86 and the RECP C augmented by additional polarization and diffuse 
functions optimized for the neutral Aus cluster [13]), 3.06 eV [18], and 3.09 eV [30]. 
Obviously, such comparison is fairly satisfactory. Let us recall that the neutral parent 
of Au5~ is the ground-state cluster Aug whereas the neutral parent of Aul^^~ is the 
higher-energy isomer Au|^^. On the neutral PES, the transition state between Aug and 
Au|^^ involves the breaking of two Au-Au bonds and its barrier is estimated to be as 
large as 15 kcal-mol~^ (see Table 2). Since a similar transition barrier exists on the 
anionic PES, Aug" is metastable and that may be why it is this isomer that has been 
observed experimentally. It would be interesting to detect experimentally the actual 
ground-state anion Au5^^~ reported here for the first time (see note [35]). 

We also demonstrate that two electrons can be simultaneously attached to Aug^^ 
to form the dianion Au^^^'^^ shown in Figure 1. Its Au-Au bond lengths are collected 
in Table 1. This dianion is stable relative to Aug^^ by 65.6 kcal-mol~^ (65.8 kcal-mol~^ 
after ZPVE) although it still remains unstable with respect to the corresponding anion 
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Aug''^''^" by 28.8 kcal-mol~^ (28.6 kcal-mol~^ with ZPVE). It is responsible for the second 
adiabatic electron affinity (see note [32]). Such dianion may therefore exist only in 
solvent or condensed phase (see Ref. [36] for current review) rather than in the gas 
phase implying that there is no contradiction with the statement that the smallest 
stable dianion is Au^^ [37]. 

Ionization of all low-energy clusters Aug"^^ yields the same well-known planar X- 
shape cation Aug"*" of symmetry [14] displayed in Figure 1. Its bond lengths, 
evaluated at the B3LYP/X (X = A, B, C) computational levels, are slightly larger, by 
0.08 A, compared to those obtained by Gilb et al. [14] using the DF BP86 in conjunction 
with the RECP C augmented by additional polarization and diffuse functions optimized 
for the neutral Aus cluster (see Table 1). The enthalpy AH/ of formation of Auf"^ from 
the ground-state pentamer cluster Aug is estimated being equal to 7.16 [7.13; 7.09] eV. 
It is consistent with the B3LYP adiabatic ionization energy lEf "^^^^(Aug) = 7.16 [7.13; 
7.09] eV. These estimates are close to the experimental data on IE^^*'*(Au5) = 8.00 
(upper hmit) [21, 22]; 7.61 ± 0.20 eV [40] and earher theoretical data: IE*''^"(Au5) = 
7.46 [6], 7.5 [18], 7.68 [21, 22], and IEf^'"-(Au5) = 7.56 (7.72) [41] and 7.78 (7.83) eV 
[9] . While the ionization pathways from Aug and Au|^ to Aug"*" are direct although the 
former is related to the AU2-AU5 bond breaking, the pathway Aug^''^ =^ Aug"*" goes via 
a transition state lying pa 15 kcal-mol~^ higher than the product and is similar to the 
one that connects Au| and Aug^^. 

The cationic PES of the five-gold clusters exhibits the first indication of the 2D - 3D 
coexistence at the B3LYP/A and B3LYP/B computational levels (see also Refs. [12b, 
41]). Together with the planar Aug"*", it contains the 3D rotamer Aug"*"' (see Figure 
1) characterized by dihedral angle ZAu2Au3Au4Aug = 82.4" and bond lengths r(Aui- 
AU2) = 2.768 and 2.629 A and r(Au2-Au3) = 2.645 and 2.629 A for the B3LYP/A and 
B3LYP/B methods, respectively. These bond lengths are very close to those in Aug"*". 
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Aug"*"' is almost isoenergetic to Aug"*". It is 0.3 kcal-mol ^ higher at the B3LYP/A and 
0.1 kcal-mol-i at the B3LYP/B after ZPVE (cf. 2.08 kcal-mor^ in Ref. [41]). The 
B3LYP/C level predicts this rotamer to be exactly energetically equivalent to Aug'*', 
though with a smaller dihedral angle, ZAU2AU3AU4AU5 — 18.0°. The difference in 
entropy at room temperature between the 2D and 3D cation of Au| varies between 
0.9 - 1.0 kcal/mol, depending on a level of computations, that together with the cor- 
responding enthalpy difference of w 0.01 eV gives a minor contribution to that of the 
B3LYP Gibbs free energies of formation of Aul,'^ and Auf"^'. 

IV. Auf 

The global-minimum energy cluster Aug is well separated from the other isomers, 
including the 3D ones in particular (see Figure 2 and Table 3). Structurally, Aug 
retains the triangular motif and is obtained from Aug by adding another triangle to 
the AU2-AU5 bond of Aug (Figure 1). It also retains the parent nonpolarity. Note that 
the BP86/A and B3LYP/B predict almost equivalent geometries of Aug (see Table 3 
for details) . The addition of one gold atom to Aug elongates its AU2-AU5 bond of by ^ 
0.13-0.16A(as the bond between 4-coordinated atoms) and shortens the bonds Au2- 
Au3 and AU4-AU5 by 0.04 - 0.05 A. Energetically, the binding energy of Aug relative to 
the asymptote Auf + Aui is equal to 60.4 (60.0 after ZPVE) [61.3 (60.9); 59.4 (59.0)] 
kcal-mol"^ (cf. 50.3 kcal-mol"^ in Ref. [9] and 63.9 (54.2) kcal-mol'^ in Ref. [29]). 

As expected, since it is a closed-shell cluster characterized by a larger binding energy 
of the paired valence electrons compared to the single valence one, Aug is more inert 
than Auf. This is confirmed by larger Ehl of 2.23 eV (BP86/A) (B3LYP: 3.48 [3.35; 
3.48] eV). The BP86/A magnitude of the Ehl is very close to the experimental value 
of 2.5 eV reported by Taylor et al. [42, 33] (compare with the earlier theoretical 
data: 2.05 eV [4, 6], 2.06 eV [9] and 1.74 eV [29]; see also the recent work [43]). The 
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IR strongest Au-Au stretching mode (doublet; 6 [5; 5] km-mol ^) is predicted at 167 
[174; 167] cm-^ 

Under ionization Aug retains its planar triangular structural form [14]. Figure 2 
displays the cation Aug"*" whose geometry is compared to that of Aug in Table 3. In- 
specting Table 3, it is worth mentioning that the B3LYP/B bond lengths of Aug"*" 
practically resemble those calculated in Ref. [14], using the DF BP86 and the RECP 
C augmented by additional polarization and diffuse functions optimized for the neu- 
tral Au3 cluster, within the margin of 0.02 - 0.06 A. The adiabatic ionization energy 
IEf3^^^(Au^) = 8.46 [8.28; 8.33] eV is about 1 eV larger than lEf ^^^^(Au^) and sat- 
isfactorily agrees with recent experimental data, IE'^^^*(Au6) = 8.80 eV (upper limit) 
[35, 36], and IE*'*^°''(Au6) = 7.60 [11], 8.55 [9], 8.17 [18], and 8.37 eV [21, 22]. 

We do not find any new clusters on the low-energy section of the PES of Aug. 
Notice that, as mentioned in Introduction, the pentagonal pyramid lies 0.44 eV higher 
Au^ (MP2/Hay-Wadt RECP) [7]. The energy difference increases to 0.65 eV [18] while 
using the MP2 method in conjunction with the 19-electron Stuttgart-type RECP C 
defined in Ref. [18] and to 0.89 eV at the BPW/B computational level. However, new 
ones are identified on the respective cationic and anionic PESs. There are two low- 
energy structures on the anionic PES (see Figure 2 and Table 3) . One of them is the 
known planar hexagonal structure Aug~ [4, 18, 33, 42, 44] structurally resembling Aug. 
Its VDE(Au^-) = 2.27 [2.18; 2.20] eV and EA„(Au^) = 2.16 [2.09; 2.09] eV (BP86/A: 
2.43 eV satisfactorily agree with the experimental data: VDE'=^^'*(Au6 ) = 2.29 eV [4], 
2.0 eV [33], 2.13(2) eV [30] (VDE*'*^''^(Au6 ) = 2.35 eV (obtained with the BP86 and 
the RECP C augmented by additional polarization and diffuse functions optimized for 
the neutral Aus cluster [13]), 2.24 eV [18], 2.29 eV [30]) and EA^^^'*(Au6) = 2.0510 ± 
0.0020 [45]; 2.016 ± 0.050 [33]; 2.06(2) eV [30] (EA*^^''^(Au6) = 2.22 eV [30]). The gap 
Ehl^ of Au^- significantly decreases to 0.4 eV (BP86/A) (B3LYP gives 1.11 [1.15; 1.10] 
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eV) compared to the neutral parent whereas the E^^j^ holds approximately at the same 
value of 2.09 eV (BP86/A) (B3LYP: 3.30 [3.27; 3.27] eV). The electronic energy and 
enthalpy of the reaction Anl,^ + Au =^ Aug^ are equal to 36.7 (36.4 after ZPVE) and 
42.6 kcakcal-mol"^ (B3LYP/A). These values are within the experimentally measured 
margin < 62.5 ± 9.2 kcal-mol"^ [46]. 

The other anion Aug^~ (Figure 2 and Table 3) is directly related to the most stable 
cluster Au5^^~ on the anionic PES of AU5. Actually, adding a single Au atom to the 
Au3 atom of Aul^^~ in Figure 1 yields Auq^~ with a gain of energy of 26.9 (26.7) 
kcal-mor^ (B3LYP/A). The anion Aug^" is however less stable, viz., by 2.5 (2.5) [5.8 
(5.7); 4.4 (4.3)] kcal-mol'^ (BP86/A: 6.7 kcal-mol'^) than Au^". Its HOMO-LUMO 
gaps amount to 1.91 eV (spin-up, BP86/A; B3LYP: 2.80 [2.82; 2.84] eV) and 0.86 eV 
(spin-down, BP86/A; B3LYP: 1.84 [1.86; 1.74] eV). Notice finally that the reaction 
pathway Aug + e~ ^ Aug^~ is indirect and occurs through a barrier. 

There are three low-energy clusters Aug"*", Aug^^"*", and Aug^"^ on the cationic PES 
shown in Figure 2 (see also Table 3 for their geometrical parameters). In fact, they are 
all cationic isomers of the neutral Aug resulting from the barrierless (Aug"^) or barrier- 
type ionization pathways. The former two are the known planar clusters [14] although 
the used B3LYP computational level shghtly favors Aug''^''^'^ over Au^+ (by 0.9 (1.0) [0.3 
(0.3); 1.4 (1.4) kcal-mol"^]; an opposite, as though negligible trend, of -0.03 kcal-mol~^, 
is predicted in Ref. [14] for the DF BP86 with the RECP C augmented by additional 
polarization and diffuse functions optimized for the neutral Aua cluster). In terms of 
the Gibbs free energy, the present order holds (with the corresponding difference of 
about 1 kcal/mol) at room temperature. The third one, Aug^"*", is 3D and somewhat 
less stable (by 5-7 kcal-mol~^) than Aug"^. Since this energy difference is close to the 
estimated energy error of the DF methods, this indicates that the 2D - 3D demarcation 
line appears at lower energies for the positively charged even-sized gold clusters than 
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for their neutral even-sized parents. 



V. Auf 

Adding a single gold atom to Aug and holding the planar triangular motif rule, 
one obtains the known global-minimum energy structure Au^ in the electronic '^A' 
state of Cs symmetry characterized by a stabihzation energy equal to 24.6 (24.4) [26.8 
(26.6); 24.6 (24.4] kcal-mol~^. This structure is shown in Figure 3 and its optimized 
geometrical parameters are gathered in Table 4 including those computed within the 
BP86/A method. Auy is an open-shell cluster, shghtly polar (0.4 [0.3; 0.3] D), having 
the following spin-up and spin-down HOMO-LUMO gaps: 1.32 and 1.04 eV (BP86/A), 
respectively (the B3LYP: 2.24 [2.19; 2.20] and 2.25 [2.20; 2.27] eV). 

Our computations show that seven is the smallest number of gold atoms leading to 
a cluster whose neutral PES contains a low-lying 3D structure. This is the novel 3D 
Auy^ cluster that is higher in energy than Auy by about 5-6 kcal-mol~^, viz., 5.6 (5.5) 
[6.5 (6.5); 5.7 (5.4)] kcal-mor^ (the BP86/A gives 5.0 (5.0) kcal-mol"^) [47]. Figure 3 
shows that the three dimensionality arises from capping one of edge triangles of Aug 
with a gold atom (called 'edge-capped pyramid' by analogy with Ref. [3]). The 3D 
structure Au^^ has twice the polarity (0.8 [0.9; 0.9] D) of Au^. Its gap Ehl^, equal to 
0.92 eV (BP86/A) (B3LYP predicts a wider gap of 1.83 [1.74; 1.78] eV), is responsible 
for its higher reactivity. The spin-down gap Ehli is larger, namely, 1.79 eV at the 
BP86/A level (cf. the B3LYP estimates it as equal to 3.04 [3.11; 3.10] eV). 

The 2D - 3D coexistence between the two structures, Auy and Auy^, appears more 
clearly on the cationic PES shown in Figure 3 (see also Tabic 4). Upon ionization, 
the 2D cluster Auy becomes the 3D cation Auy^"*" via the 2D barrier Auy"*"*^, which 
is structurally similar to its parent neutral. Auy^"*" is similar to the 3D neutral Au^^ 
(Figure 3). The transition-state linker Auy'*'*'' lies 4.8 kcal-mol~^ higher (B3LYP/A). 
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This is one remarkable facet of the cation Auy^"*". The other is that ionization of Au^^ 
directly yields Auy^"*". Such a 2D - 3D geometry change upon ionization is reported 
in this work for the first time. The Auy =^ Auy^** =^ Auy''^^ pathway determines the 
ionization energy lEf 3-^^^(Auf) = 6.99 [6.94; 6.91] eV which is smaller than lEf ^^(Aur) 
= 7.80 eV [35, 36] and rather close to the theoretical value IE*'*^'"'(Au7) = 7.18 given 
inRef. [18]. 

Upon electron attachment, the low-energy 3D cluster Au^^ found in this work ex- 
hibits a 3D to 2D conversion that leads to the novel planar structure Auy^". It is 
shown in Figure 3 (see also Table 4 for its geometry) and characterized by rather large 
VDE(Au7^") = 4.50 [4.34; 4.42] eV. The latter clearly indicates that a 3D structure for 
the anion is less stable than the less strained 2D one. The HOMO-LUMO gap of Auy^" 
is equal to 1.72 eV (BP86/A) (the B3LYP: 2.85 [2.77; 2.79] eV). However, Auf" does 
not occupy the global minimum on the anionic PES of Auy, rather it lies above the 
well-known ground-state tricapped square cluster Auy" [3, 30] by 4.6 (4.5) [7.9 (7.7); 
6.0 (5.9)] kcal-mol-^ (the BP86/A estimate is 10.5 (10.3) kcal-mol"^). Due to a large 
entropy effect, the difference in the B3LYP Gibbs free energies between Auy^" and 
Auy", evaluated at T = 298 K, is rather small and equal to 1.9 [4.7; 2.2] kcal-mol~^, 
reflecting the floppy nature of the bottom of the anionic PES. 

The anionic PES exhibits two energetically almost degenerate 2D structures: Auy" 
and Au7^^~ shown in Figure 3 (see Table 4 and note [48]). The latter is a novel structure 
obtained by adding a single gold atom to one of the edge atoms of Aug". Actually, it 
lies 2.5 (2.5) [1.7 (1.7); 2.4 (2.4)] kcal-mol'^ (the BP86/A: 0.9 (0.9) kcal-mor^) below 
Au^~. Its VDE(Au^^^-) = 3.87 [3.74; 3.74] eV and its HOMO-LUMO gap is equal to 
1.35 (BP86/A) (whereas the B3LYP gives 2.17 [2.17; 2.12] eV). 

The parent neutral cluster of Au^''^''^^ is the honeycomb Au^''^''^ in the state ^i?2g of 
symmetry [10]. On the neutral PES, Au^^^ lies however 9.0 (8.8) [7.7 (7.5); 9.0 (8.8)] 
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kcal-mol ^ above the ground state (the BP86/A value is 8.1 (7.9) kcal-mol ^; cf. with 
7.0 kcal-mol~^ in Ref. [10]; see also note [49]). In contrast to Au^, Au^^^ is nonpolar 
and its spin-split HOMO-LUMO gaps are equal to 0.18 and 1.73 eV (BP86/A) (cf. 
with the B3LYP predictions: 1.01 [0.95; 0.97] and 2.76 [2.76; 2.78] eV and the value of 
0.87 eV reported in Ref. [29]). 

Compared to Au^y^~, Auy" is more reactive due to a smaller Ehl gap = 1.18 eV 
(BP86/A) (the B3LYP: 1.95 [1.94; 1.89] eV) and is characterized by an EA„ = 3.40 
(3.40) [3.29 (3.29); 3.28 (3.28)] eV (3.62 eV for the BP86/A) and a VDE = 3.50 [3.39; 
3.38] eV. Both these computed values are fairly close to EA^^p*(Au7) = 3.415 ± 0.050 
[33] and 3.40(3) eV [30] and VDE^^P*(Auf) = 3.46(2) eV [3, 30] and 3.38 eV [33], 
respectively. They also fairly agree with previous theoretical data: EA*'^^°^(Au7) = 
3.43 eV [30] and VDE*'*'=''^(Auf ) = 3.55 eV (obtained with the BP86 and the RECP C 
augmented by additional polarization and diffuse functions optimized for the neutral 
Au3 cluster [13]), 3.43 eV [18], and 3.46 eV [30]. 

Let us return to the cationic PES of Auy sketched in Figure 3. Two novel different 
structures are found on its lowest-energy portion. One, as described above, is the cation 
Auy"*" and the other is a 3D honey comb- type Auy^^"*" (see Table 3 and note [50]). The 
energy difference between Au^^^+ and Au^+ is 0.1 (0.2) [2.9 (2.9); -1.9 (-2.0)] kcal-mol"^ 
[51] meaning that the computational level B3LYP/C slightly favors the cation Auy^^^. 
As shown in Figure 3, this cation Auy^^"*" is not perfectly planar at the B3LYP/A and 
B3LYP/B computational levels which predict the dihedral angle ZAU2AU3AU7AU6 = 
31.6° and 31.5°, respectively (B3LYP/C gives only 0.2° and BP86/B does 27.5°). In 
summary, it is found for the first time that there may be no 2D clusters in the low- 
energy section (likely up to the energies of ca. 15 kcal-mol~^; see Ref. [10] for the 
corresponding neutrals) of the cationic PES of seven gold atoms. 

To end this Section, let us mention another neutral structure Auy^ in the ^Ai state 
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of symmetry (Figure 3 and Table 3) lying 8.5 (8.6) kcal-mol~^ (B3LYP/A) above 
Auy^^ (cf. with 8.3 kcal-mol~^ reported in Ref. [10]). It possesses spin-up and spin- 
down HOMO-LUMO gaps of 1.52 and 1.50 eV, respectively. The related transition- 
state double-rhombus structure Auy^** is also shown in this Figure 3 and determines a 
rather high barrier of falT.T kcal-mol~^ governing the breaking of the AU2-AU5 bond. 

VI. Auf 

It is known that the Aug cluster is the smallest catalytically active size of gold 
clusters [le]. 

A binding of a single gold atom Aui to the AU4-AU5 bond of the Auy cluster leads 
to the planar and nonpolar tetracapped square cluster Au| (Figure 4) which is the 
global minimum on the neutral PES of eight atom gold clusters. Such a structure has 
recently been reported by Bonacic-Koutecky et al. [10, 21] as the true global minimum, 
assigned to the ^Aig state of symmetry, in contrast to earlier findings. The binding 
energy of Aug, taken relative to the infinitely separated Auy and Aui, is equal to 56.0 
(55.6) [56.8 (56.5); 54.6 (54.2)] eV. The Ehl gap width amounts to 1.31 eV (BP86/A) 
(the B3LYP gives 2.77 [2.72; 2.67] eV). 

Structurally, Aug is composed of two kinds of coordinated gold atoms: the four 
edged atoms are 2-coordinated while the other four (central square) are 4-coordinated. 
The Au-Au bond length between 2- and 4-coordinated atoms is about 0.15 - 0.17 A (see 
Table 4) shorter than between each pair of 4-coordinated atoms. Each pair of opposite 
atoms belonging to the central square are unbounded, forming thus a sort of hole since 
the corresponding distances of ca. 4.0 A exceed twice the van der Waals radius of gold 
(3.32 A). It is worth noticing that the formation of Aug from Auy -|- Aui requires to 
break an Au-Au bond of length f« 3.0 A (< 3.32 A) in the internal rhombus of Auy. 

The two following low-energy structures, Aug^ and Aug^^, have almost the same 

17 



energy but differ by their 2D - 3D character (see Figure 4 and Table 2). They he about 
11 - 12 kcal-mol-^ above Au^ : AE(Auf - Au^) = 12.5 (12.4) [10.8 (10.7); 11.0 (10.9)] 
kcal-mol-i (and 9.7 (9.6) kcal-mol-i for the BP86/A) and AE(Au^^^ - Au^) = 12.7 
(12.6) [13.9 (13.8); 12.8 (12.7)] kcal-mor^ and are of the same polarity («i 0.9 D). One 
of them, Aug^, is a planar C2v structure that possesses a sixfold coordinated atom at 
the center. The other, Au|^^, is a novel 3D cluster formed by a Au2 dimer bonded to 
the six-gold most stable cluster Aug along its symmetry axis. It possesses four fivefold 
coordinated atoms of gold. The difference between the two clusters is also manifested 
in the width of ^hl'- while for the 3D cluster Au|^^, it is closer to that of Aug, viz., 
1.86 eV (BP86/A) (B3LYP: 2.95 [2.85; 2.94] eV), its 2D counterpartner has a narrower 
one, i. e., 0.87 eV (BP86/A) (B3LYP: 1.76 [1.71; 1.72] eV). These two structures are 
shghtly below (ca. 3-4 kcal-mol~^ depending on the computational level) a pair of less 
stable novel 3D clusters, Aug^ and Au^ shown in Figure 4. Aug^ is formed by binding 
a gold dimer to Aug along its 3-atom bond whereas in Au^ a gold trimer parallel is 
attached parallel to Aug (see Figure 1). 

While the 3D cluster Aug^ is rather far from the global minimum of the neutral 
PES, on the cationic one, Aug^"*" competes with the 2D Aug"*" and Aug^"*", and the three 
isomers he at the very bottom of the PES (see Figure 4). More precisely, Aug^"*" is only 
0.1 (0.1) [2.3 (2.3); 1.2 (1.2)] kcal-mor^ above Au^+ which in turn is only 1.7 (1.7) [0.6 
(0.6); 0.4 (0.4)] kcal-mol"^ lower than Au8^+ [52]. This is a manifestation of the 2D 

3D transition that occurs on the cationic PES of eight gold atoms. The optimized 
geometries and some other properties of these three cationic clusters are presented in 
Table 5. Notice that the B3LYP/B bond lengths of Aug^"*" are in a good agreement 
(within the margin of 0.03 - 0.04 A) with the bond lengths obtained in Ref. [14] usmg 
the BP86 and the RECP C augmented by additional polarization and diffuse functions 
optimized for the neutral Aus cluster. The lEa(Aug) = 7.88 [7.71; 7.72] eV correlates 



18 



with the experimental upper hmit of 8.65 eV [35]. The IE„(Au^^) = 7.41 [7.27; 7.26] 
eV. 

That Aug^ is the most stable anionic cluster is a well established fact [13]. Its 
structure is displayed in Figure 4. Two other low-energy clusters: the planar Aug''^" 
which is 7.1 (7.0) [6.9 (6.8); 6.6 (6.6)] kcal-mol"^ higher than Aug" and the three- 
dimensional Aug^^" lying 21.1 (20.9) [23.4 (23.2); 21.0 (20.8)] kcal-mol"^ higher Aug" 
are shown in Figure 4 as well (see also Table 5). The HOMO-LUMO gaps of the open- 
shell cluster Au^- amounts to 1.01 eV (BP86/A) (cf. B3LYP: 1.69 [1.73; 1.68] eV) 
(spin-t) and 1.31 eV (BP86/A) (cf. B3LYP: 2.51 [2.53; 2.44] eV) (spin-j), respectively. 
The computed adiabatic electron affinity of Aug" is equal to 2.78 [2.74; 2.72] eV in 
agreement with the experimental magnitude of 2.764 ± 0.050 eV [33] and 2.73(2) eV 
[30] and the theoretical one of 2.88 eV [30]. The VDE of Aug" amounts to 2.88 [2.82; 
2.82] eV which is also close to VDE^^^'*(Aug") = 2.79 [33] and 2.79(2) eV [30] and 
YQgtw^A^-) ^ 3 eV (obtained with the BP86 and the RECP C augmented by 
additional polarization and diffuse functions optimized for the neutral Aua cluster [13]) 
and 2.94 eV [30]. Note that the 3D anion Au|^^~ is stable under attaching e~, that 
is, the dianion Aug^^^~, which becomes 2D and whose structure is displayed in Figure 
4 (see also Table 5), is energetically lower than Aug''^^" by 14.5 [11.8; 12.2] kcal-mol~^ 
and therefore exists as a long-lived metastable dianion. 

VII. The Low-Energy Section of the PES of Aug 

As we have seen in the previous Sections, the planar triangular motif works only 
for the neutral and anionic ground-state clusters Aug_g and not for the cationic ones. 
Does it mean that it is applicable for all larger clusters or is there some threshold size 
of gold clusters for which this structural rule is no longer valid? 

We show in this Section that n = 9 is in fact such a threshold. The structures of the 
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low-energy nine gold atom clusters, displayed in Figure 5 (see also Table 6), supports 
this conjecture: three of them, a bicapped hexagonal ring Aug, Aug^^, and Aug^ are 
2D whereas Aug^ is 3D. The clusters Aug, Aug-'^ and Aug^ are novel and found in the 
present work for the first time (cf. the clusters in Refs. [4, 10]; see notes [53, 54]). 
The reported four structures of Aug can be considered as nearly isoenergetic due to 
the definition of the accepted of the DF/RECP methods for the gold clusters (see Ref. 
[10] and Section II). 

Comparing the neutral PESs of 6 - 9 gold atoms, we conclude that the three pla- 
nar low-energy structures Aug, Aug'^''^ and Aug^ originates from Aug, Auy or Aug by 
adding one, two or three planar triangles, respectively. On the other hand, the almost 
isoenergetic 3D structure Aug^ is directly obtained from Aug by adding (capping) one 
gold atom to its central square and forming thus four new Au-Au bonds, instead of two 
as demanded by the planar triangular motif rule. This means that n = 9 is actually 
the minimal number of gold atoms for which the above rule is no longer fully valid. 
Therefore, the neutral PES of nine gold atoms markedly demonstrates the 2D - 3D 
coexistence on its low energy section. 

VIII. Summary and Outlook 

A variety of novel low-energy 2D and 3D structures have been determined in the 
present work on the neutral, cationic and anionic PESs of small gold clusters Au5<„<8 
and on the neutral PES of Aug. They provide new insights on the 2D - 3D transition 
in gold clusters, particularly on their well-known odd- vs. even-sized patterns. The 
neutral PES of Aus shows the presence of two new planar forms, Aug^^ and Aug^, 
fiUing thus the gap of « 21 kcal-mol~^ between the 2D ground-state W-form and the 
3D pyramid discussed in Ref. [10]. On the PES of Au7, this gap is significantly reduced 
to ~ 6 kcal-mol"^ due to the existence of the novel 3D cluster Au^^. For n = 9, the gap 
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is essentially closed: the novel 3D structure Aug''^ coexists on almost equal energetical 
footing with three 2D structures, two of which, including the slightly most stable one, 
Aug, are new as well. 

As far as the 2D - 3D coexistence is concerned, the small even-numbered closed- 
shell clusters (n < 10) behave quite differently from the small neutral odd-numbered 
open-shell ones. This is particularly true for the PES of Aug where the lowest-energy 
3D cluster Aug^^ is separated from the 2D ground-state one by ca. 13 kcal-mol~^. We 
suggest that this is due to the stronger s — d hybridization in Au2n than in Au2n+i, 
at least for 2 < n < 4. Another feature related to the 2D - 3D transition in small 
odd-numbered clusters of gold emerges from the present work. This is the behavior 
of the gap EffL- The spin-up gap of the ground-stable 2D clusters narrows with 
increasing n from 1.57 eV (Au^) to 1.32 (Au^) and 0.90 eV (Au^; the Ehl^ of Au^^ = 
0.22 eV). On the contrary, for n = 7 and 9, the gap of the lowest-energy 3D clusters 
Au^^ and Au^^ widens from 0.92 to 1.32 eV. Note that E^^^^ of the 3D cluster Au^^ 
is smaller by 0.4 eV than that of the 2D one Auy. This implies that the former has 
a higher catalytic activity despite its smaller binding energy and therefore deserves to 
be experimentally tested. The spin imbalance of the HOMO-LUMO gap, Eh^ - Ehli, 
decreases from 0.57 eV at n = 5 to 0.28 eV at n = 7 and 0.14 eV at n = 9 (Aug). 

Novel neutral clusters also lead to interesting features on the charged PESs. One 
example: the neutral cluster Au|^^ gives rise on the anionic PES to the anion Aug^^" 
which has been found to be the most stable one, in contrast with previous experimental 
and theoretical works reporting that the ground anionic state is the W-shape structure 
Aug". This discrepancy can be explained by the presence of a pair of similar transition 
states, one between the ground-state W-form Au| and higher-energy cluster Au|^^ on 
the neutral PES and another between Aug" and Aug^^" on the anionic PES. These 
transition states govern the breaking of two Au-Au bonds and are therefore charac- 
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terized by rather high barriers (of at least 15 kcal-mol~^ on the neutral PES) which 
results in a lifetime of Auf~ long enough to be experimentally observed. On the other 
hand, the Aug^^" isomer, though less reactive than Au5~ due to a larger Fj^l, would be 
interesting to detect as well for another two reasons. First, as demonstrated in Section 
III, it gives rise to the long-hved metastable 5-gold dianion whose structural pattern is 
still present in the 8-gold dianion Aug^^^". Second, the bonding of a single gold atom 
to the atom Aus of Au|^^~ results in the formation of the anionic cluster Aug^~ which 
hes only 3-5 kcal-mol~^ higher the ground-state anion Auq~. Both of them are built, 
by analogy with the triangular motif used for neutrals, via the bonding of Au to the 
respective 5-gold anions. 

As established above, for neutral clusters, the 2D - 3D coexistence takes place only 
for the odd-sized open-shell gold cluster Aug. However, it already develops for the 
cationic odd-sized closed shell clusters Au^ and Auy , although the ways it does for 
these two n are essentially different. For n — 5, the 2D - 3D transition emerges from 
the existence of two almost isoenergetic rotamers: the 2D Aug"*" and the 3D Aug^' 
which lacks its neutral 3D parent. On the other hand, the case n = 7 is directly related 
to the ionization of the ground-state neutral 2D cluster Au^. As noted in Section 
V, the ionization weakens the s — d hybridization and thus makes the 2D triangular 
motif rule invalid. We suggest that this process occurs via the intermediate planar 
transition-state hnker Auy"*"*^ which relaxes into the 3D cluster Auy^"*". In contrast to 
the 3D cation Auf"*"', Auy^"*" possesses its 3D parent image Au^^ on the neutral PES 
which however, as already discussed above, is less stable than Au^. Interestingly, the 
3D Zn-analogue of Auy^"*", where the Zn atom replaces Au4 (or Au5, see Figure 3), is 
only 3.2 kcal-mol~^ (B3LYP/A) higher the most planar stable cation AueZn^ reported 
in Ref. [20a] (the cluster [6a+] therein), that is, it falls within the DF/RECP error 
margin. 
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In addition, the present computational study clearly demonstrates that another 
cationic cluster Au^^^'^ (see Figure 3), originating from the neutral planar honeycomb 
structure Auy^^, is also 3D (at the B3LYP/A and B3LYP/B computational levels). As 
reported in Section V, it is almost isoenergetic to Auy^"*" within the DF/RECP error 
margin. Moreover, the cation Auy^"*", separated from Auy^"*" by about 9.6 kcal-mol~^ 
(B3LYP/A), is 3D too (see Figure 3). Altogether, this means that in the low-energy 
section of the PES of Auy" , 3D structures are favored over the 2D ones up to energies 
of ca. 15 kcal-mol"^ (see the neutrals in Ref. [10]). In summary, the low-energy portion 
(within of ca. 15 kcal-mol"-*^ from the bottom) of the cationic PES of seven gold atoms 
is characterized by a complete absence of 2D clusters. Strickly speaking, the concept 
of the 2D - 3D coexistence might then be only vahd beyond these energies. 

The 2D - 3D coexistence in even-sized open-shell cationic clusters of gold is quite 
different from that of the odd-sized closed-shell cationic clusters and that of the even- 
sized closed-shell neutrals as well. As an example, the novel 3D cation Aug^"*" lies 
only 5-7 kcal-mol^^ higher than the almost degenerate ground-state cations Aug^ and 
Au^^^+. At n = 8, the energy offset between the 2D ground-state cluster Aug^"*" and 3D 
one Aug^"*" reduces to f« 1 kcal-mol~^ which falls within the DF/RECP computational 
error. One may therefore anticipate that the even n — 10 might show the clearcut 2D 
- 3D coexistence at the bottom of the PES. 

The final comment is about the anionic PESs of gold clusters. It has been demon- 
strated in the present work that the studied anions hardly manifest the 2D - 3D co- 
existence. This conclusion is in agreement with the early findings by Landman and 
co-workers [4, 30]. The most stable 3D cluster which has been reported in this work for 
the first time is Aug^^". However, it is strongly (by about 21 - 23 kcal-mol~^) separated 
from the 2D ground-state cluster Aug~. 
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Figure Legends 



Figure 1. The low-energy sections of the neutral, cationic and anionic PESs of Au5 
and the structure of the long-hved dianion Aug^^^". 

Figure 2. The low-energy portions of the neutral, cationic and anionic PESs of six 
gold atoms. 

Figure 3. The low-energy portions of the neutral, cationic and anionic PESs of 
seven gold atoms. 

Figure 4. The low-energy portions of the neutral, cationic and anionic PESs of Aug 
and the structure of the long-hved dianion Aug^^^". 

Figure 5. The low-energy portion of the neutral PES of nine atoms of gold. The 
bond lengths of the displayed clusters are listed in Table 6. 
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Table 1: The bond lengths (in A) of the low-energy neutral and charged clusters of five gold atoms. 





B:!LYr/A 


B:!LYr/D 


D:iLYP,'r 


BP.-:ii,'A 


rAV')l/A 


Dill 11 


r(Aui-Au2) 
r(Aui-Au3) 
r(Au2-Au3) 
r(Au2-Au5) 


2.872 
2.698 

2. no 

2.777 


2.831 
2.678 
2.722 
2.740 


2.863 
2.702 
2.762 
2.777 


2.837 
2.667 
2.723 
2.734 


2.829 
2.666 
2.721 
2.730 


2.78 [4]; 2.81 [7]; 2.704 [8] 
2.64 [4], 2.72 [7], 2.566 [8] 
2.75 [6]; 2.614 [8] 
2.79 [Dj; J.DU9 [esj 


Au^-(lAiC2^,) 
r(Aui-Au2) 
r(Aui-Au3 ) 
r(Au2-Au3) 
r(Au2-Au5) 


3.001 
2.679 
2.881 
2.676 


2.906 
2.666 
2.831 
2.831 


2.959 
2.685 
2.877 
2.877 






2.82 [4]; 2.725 [8] 
2.65 [4]; 2.566 [8] 

2.672 [8] 

2.582 [8] 


AuJ+(^B2. D-^r,) 
r(Aui-Au2 ) 
r(Au2-Au5) 


2.765 
2.643 


2.736 
2.628 


2.767 
2.650 






2.684 [14] 
2.588 [14] 


5 

r(Aui-Au2 ) 
r(Au2-Au3) 


2.768 
2.645 


2.629 
2.629 










Au^^(2B2^ D2h) 

r(Aui-Au2) 

r(Au2-Au5) 


2.714 
2.751 


2.689 
2.722 


2.716 
2.751 


2.684 
2.723 




2.609 [8] 


r(Aui-Au2 ) 
r(Aui-Au3) 
r(Aui-Au4) 
r(Au4-Au5) 


2.619 
2.731 
2.793 
2.619 


2.603 
2.707 
2.742 
2.601 


2.631 
2.729 
2.793 
2.629 


2.594 
2.684 
2.742 
2.586 






A III- 

5 

r(Aui-Au2) 
r(Aui-Au3) 
r(Aui-Au4) 


2.628 
2.719 
2.922 


2.613 
2.697 
2.697 


2.635 
2.718 
2.718 








A III2.- 

Aug 

r(Aui-Au2) 
r(Aui-Au3) 
r(Aui-Au4) 


2.740 
2.867 
2.734 


2.707 
2.818 
2.715 


2.728 
2.848 
2.746 








r(Aui-Au3) 
r(Aui-Au4) 

r(Au2-Au5) 
r(Au3-Au4) 
r(Au3-Au5) 


2.872 
2.643 
2.650 
2.744 
2.643 


2.814 
2.640 
2.626 
2.701 
2.617 


2.839 
2.645 
2.662 
2.769 
2.657 









Table 2: The relative energies (in kcal-mol^^) AE and ZPVE-corrected ones AEzpve 
(in parentheses) and the HOMO-LUMO gaps Ejjl (in eV) of the neutral, cationic and 
anionic gold clusters Au5<„<9. Notice that En^i-^'^i^) — ^-^^ ®^ [29]- 





D:iLYr,'A 


AE (AE 
D:!LYr , D 


D:!LYr (■ 


Eirsfi, A 


D:iLYP;'A 


D:!LYr , D 


^HL 1 

D:!LYr , (■ 


Dr<s(i/A 


AU5 

Au'r 

Aut' 

aJ" 
aJv 


0.0 
7.28 (7.27) 
13.88 (13.80) 
18.56 (18.40) 


0.0 
8.17 (8.15) 
15.93 (15.83) 
20.51 (20.34) 


0.0 
7.99 (7.97) 
14.92 (14.81) 
19.37 (19.18) 


0.0 

8.04 (8.03) 
15.16 (15.15) 


2.49 2.17 
1.44 2.84 
2.68 1.38 
1.16 2.36 


2.46 2.17 
1.50 2.73 
2.65 1.42 
1.12 2.46 


2.47 2.17 
1.43 2.98 
2.62 1.34 
1.24 2.10 


1.58 1.01 
0.77 1.31 
1.83 0.50 


Au 


0.0 
7.21 (7.19) 


0.0 
5.82 (5.84) 


0.0 
6.30 (6.30) 


0.0 


2.88 
1.21 


2.39 
0.93 


2.41 
1.09 




AU6 

Aug 










3.48 


3.35 


3.45 


2.23 


Au- 
Au„ 
Au?- 


0.0 
2.62 (2.46) 


0.0 
5.78 (5.66) 


0.0 
4.36 (4.25) 


0.0 


1.11 3.30 
2.80 1.84 


1.15 3.27 
2.82 1.86 


1.10 3.27 
2.84 1.74 


0.40 2.09 
1.91 0.86 


Au7 

Au'. 

aJ' 

Aul" 


0.0 
6.66 (5.51) 
8.99 (8.82) 


0.0 
6.46 (6.41) 
7.69 (7.52) 


0.0 
5.73 (5.43) 
8.15 (7.98) 


0.0 
5.08 (5.02) 
8.07 (7.87) 


2.24 2.25 
1.83 3.04 


2.19 2.20 
1.74 3.11 


2.20 2.27 
1.78 3.10 


1.32 1.04 
0.92 1.79 


Au" 

Aui!- 

Au^- 

, III- 
Au^ 


0.0 
4.56 (4.45) 
-2.51 (-2.50) 


0.0 
7.87 (7.70) 
-1.71 (-1.73) 


0.0 
6.01 (5.S7) 
-2.44 (-2.44) 


0.0 
2.85 
-0.87 (-0.87) 


1 . 95 
2.77 
2.17 


1.94 
2.79 
2.17 


1.89 
1.72 
2.12 


I.IS 
1.35 


Aug 

Aul 
Au?' 


0.0 

12.53 (12.44) 


0.0 

10.77 (10.72) 


0.0 

11.01 (10.94) 


0.0 
9.67 (9.62) 


2.77 
1.76 


2.72 
1.71 


2.67 
1.72 


1.74 
0.87 


Aug 

Auq 
Au ' 
Au" 

Aur 


0.0 
1.98 (1.99) 
3.22 (3.21) 
4.36 (4.26) 


0.0 

4.55 (4.53) 
3.24 (3.03) 

4.56 (4.52) 


0.0 
4.33 (4.30) 
3.03 (3.02) 
3.85 (3.05) 


0.0 
1.72 (1.71) 
2.56 (2.56) 
3.55 (3.48) 


1.80 1.66 
2.21 1.92 
2.12 1.03 
0.91 2.39 


1.74 1.65 
2.19 1.91 
2.05 1.04 
0.92 2.33 


1.77 1.62 
2.19 1.85 
2.07 1.00 
0.89 2.40 


0.90 0.76 
1.32 1.08 
1.22 0.19 
0.22 1.35 
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Table 3: The bond lengths (in A) of Au^, Au^+, Au^", Au^^^, Au^^^+, and Au^^+. 





B3LYP/A 


B3LYP /B 


B3LYP /C 


BP86/A 


Data 


Aug( A-^ D^h) 












r(Aui-Au2) 


2.900 


2.856 


2.890 


2.856 


2.81 [4]; 2.84 [7] 


r{Aui-Au3) 


2.900 


2.856 


2.890 


2.855 


2.81 [4] ; 2.84 [7] 


r(Aui-Au4) 


2.704 


2.683 


2.706 


2.677 


2.66 [4]; 2.72 [7] 


r{Aui-Au5) 


2.704 


2.683 


2.706 


2.677 


tAl. "70 ["71 

J. 00 [4J , A. I A yi\ 


r{Au3-Au4) 


2.704 


2.683 


2.706 


2.677 


2.66 [4]; 2.72 [7] 














r(Aui-Au2) 


2.997 


2.930 


3.076 




2.915 [14] 


r(Aui-Au3) 


2.732 


2.708 


2.804 




2.752 [14] 


r(Au]^ -AU4) 




2.688 


2.694 






r(Aui-Au5) 


2.700 


2.674 


2.733 






r( Au3- AU4) 


2.704 


2.679 


2.727 


















r{Aui-Au2) 


2.810 


2.781 


2.808 


2.784 


2.76 [4] 


r{Aui-Au3) 


2.810 


2.781 


2.808 


2.784 


2.76 [4] 


r(Aui-Au4) 


2.781 


2.746 


2.778 


2.743 


2.71 [4] 


r{Aui-Au5) 


2.781 


2.746 


2.778 


2.743 


2.71 [4] 


r{Au3-Au4) 


2.781 


2.746 


2.778 


2.743 


2.71 [4] 














r(Aui-Au2) 


2.764 


2.728 


2.765 


2.715 




r(Aui-Au4) 


2.646 


2.623 


2.654 


2.614 




a4"+ 












r{Aui-Au2) 


2.700 


2.680 


2.704 






r(Aui-Au3) 


2.736 


2.703 


2.739 






r{Au2-Au3) 


2.803 


2.778 


2.802 






r{Au2-Au4) 


2.824 


2.777 


2.820 






r(Au3-Au4) 


2.868 


2.821 


2.865 


















r(Aui-Au2) 


2.720 


2.720 


2.729 






r{Aui-Au3) 


2.744 


2.726 


2.753 






r(Aui-Au4) 


2.948 


2.888 


2.946 






r{Aui-Au5) 


2.999 


2.892 


2.957 






r(Au2-Au3) 


3.079 


2.948 


2.997 






r{Au2-Au5) 


2.822 


2.783 


2.831 






r{Au3-Au4) 


2.783 


2.773 


2.794 






r(Au4-Au6) 


2.712 


2.695 


2.715 






r(Au5-Au6) 


2.722 


2.699 


2.726 
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Table 4: The bond lengths (in A) of the clusters Au^ and their charged cousins. 
The superscript a indicates the bond length equal to 2.714 A in Ref. [4]. 



B3LYP/A B3LYP/B B3LYP/C BP86/A 



B3LYP/A B3LYP/B B3LYP/C BP86/A 



A' Cs) 

AU2) 
AU4) 
Aug) 
AU4) 
AU4) 
Aue) 
Aut) 
AU5) 
Aug) 
Aug) 
Au7) 



AU2) 
AU5) 
Aug) 
AU4) 
Aus) 
AU4) 
Aus) 
Aug) 
AU7) 
Aug) 
Au7) 
Au7) 



r{Aui 
r{Aui 
r{Au2 
r(Au2 
r{Au3 
r(Au3 
r{Au3 
r{Au4 
r{Au4 
r(Au5 
r{Au6- 

Au^^ 

r{Aui- 
r{Aui- 
r{Au2- 
r{Au2- 
r(Au2- 
r(Au3- 
r(Au4- 
r(Au4- 
r(Au4- 
r{Au5- 
r(Au5- 
r(Au6- 

a4" 

r{Aui-Au2) 
r{Aui-Au4) 
r{Au3-Au7) 
r{Au4-Au5) 
r(Au5-Au6) 
r(Aug-Au7) 

a4"+ 

r(Aui-Au2)" 
r(Aui-Au4) 

Au^-(lAi C2„) 

r{Aui-Au2) 

r(Aui-Au7) 

r{Au2-Au3) 

r{Au3-Au4) 

r(Au3-Au5) 

r(Au5-Au7) 

r{Aui-Au4) 
r{Aui-Au5) 
r{Au2-Au3) 
r(Au2-Aug) 



769 


2.734 


2.763 


2.745 


756 


2.725 


2.762 


2.711 


717 


2.693 


2.718 


2.683 


780 


2.753 


2.779 


2.753 


102 


3.003 


3.065 


2.980 


868 


2.829 


2.860 


2.843 


724 


2.702 


2.728 


2.698 


712 


2.689 


2.715 


2.687 


845 


2.802 


2.840 


2.794 


695 


2.674 


2.696 


2.669 


713 


2.688 


2.713 


2.683 



706 


2.685 


2.708 


2.679 


701 


2.678 


2.703 


2.674 


706 


2.685 


2.708 


2.679 


899 


2.855 


2.890 


2.856 


899 


2.855 


2.890 


2.856 


701 


2.678 


2.703 


2.674 


966 


2.917 


2.948 


2.914 


756 


2.786 


2.816 


2.777 


901 


2.786 


2.816 


2.777 


901 


2.786 


2.816 


2.777 


756 


2.786 


2.816 


2.777 


916 


2.869 


2.923 


2.885 



840 


2.725 


2.753 


2.732 


750 


2.732 


2.805 


2.717 


910 


2.791 


2.766 


2.854 


842 


2.860 


2.860 


2.802 


756 


2.810 


2.905 


2.861 


755 


2.721 


2.751 


2.733 



2.807 
2.779 



2.773 
2.747 



2.790 
2.790 



837 


2.788 


2.830 


648 


2.640 


2.655 


730 


2.712 


2.733 


719 


2.697 


2.719 


874 


2.807 


2.862 


797 


2.753 


2.791 


692 


2.675 


2.694 


882 


2.822 


2.868 


727 


2.700 


2.729 


622 


2.608 


2.629 



2.777 
2.750 



2.783 
2.634 
2.712 
2.694 
2.802 
2.752 



2.670 
2.810 
2.695 
2.597 



Au^+*<' 

r(Aui-Au2) 
r(Aui-Aui) 
r(Au2-Au3) 
r(Au2-Au4) 
r(Au3-Au4) 
r(Au3-Aug) 
r(Au3-Au7) 
r(Au4-Au5) 
r(Au4-Au6) 
r(Au5-Aug) 
r(Aug-Au7) 

Au^^+ 

r(Aui-Au2) 
r(Aui-Au5) 
r(Au2-Au3) 
r(Au2-Au4) 
r(Au2-Au5) 
r(Au3-Au4) 
r(Au4-Au5) 
r(Au4-Aug) 
r(Au4-Au7) 
r(Au5-Aug) 
r(Au5-Au7) 
r(Au6-Au7) 

A^^ 

r(Aui-Au2) 
r(Aui-Au3) 
r(Au2-Au3) 
r(Au2-Au5) 
r(Au2-Aug) 
r(Au3-Aug) 

Au^^*» 

r(Aui-Au2) 

r(Au2-Au6) 

A III- 
Au^ 

r(Aui-Au2) 
r(Aui-Au3) 
r(Au3-Au4) 
r(Au3-Au5) 
r(Au5-Au6) 



2.710 
2.713 
2.810 
2.858 
2.903 
2.918 
2.667 
2.701 
2.782 
2.700 
2.754 



725 


2.699 


2.726 


671 


2.654 


2.675 


725 


2.699 


2.726 


876 


2.841 


2.871 


871 


2.841 


2.871 


671 


2.654 


2.675 


995 


2.934 


2.979 


861 


2.822 


2.860 


861 


2.822 


2.860 


861 


2.822 


2.860 


861 


2.822 


2.860 


861 


2.670 


2.685 



2.881 
2.795 
2.818 
2.697 
2.797 
2.695 



2.751 
2.744 



655 


2.635 


2.662 


2.625 


757 


2.725 


2.758 


2.718 


766 


2.762 


2.766 


2.763 


789 


2.762 


2.786 


2.748 


692 


2.677 


2.696 


2.672 
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Table 5: The bond lengths (in A) of the neutral, cationic and anionic clusters of Aug. 





E..>L\ r A 


B.iL\ r , E. 




Br s(), . 












r(Aui-Au2) 


2.697 


2.675 


2.698 


2.670 


r(Aui-Au5) 


2.697 


2.675 


2.698 


2.670 


r{ Au2- AU4) 


2.870 


2.826 


2.860 


2.819 


r(Au4-Au5) 


4.066 


4.003 




3.980 


Au^+{2AigD4fc) 










r(Aui-Au2) 


2.722 


2.695 


2.720 




r{Aui-Au5) 


2.704 


2.679 


2.711 




r(Au2-Au4) 


2.844 


2.805 


2.840 




r(Au4-Au5) 


3.004 


2.919 


2.990 




AuJ-(2Aij,D4h) 










r{Aui-Au2) 


2.754 


2.722 


2.752 


2.721 


r{Aui-Au5) 


2.754 


2.722 


2.752 


2.721 


r{Au2-Au4) 


2.782 


2.752 


2.780 


2.761 


r(Au4-Au5) 


3.940 


3.893 


3.931 


3.895 


Au" 

s 










r(Aui-Au2) 


2.946 


2.866 


2.924 


2.866 


r{Au2-Au3) 


2.771 


2.743 


2.766 


2.764 


r(Au3-Au4) 


2.809 


2.766 


2.805 


2.758 


r(Au3-Aug) 


2.877 


2.827 


2.866 


2.816 


r(Au5-Aue) 


2.733 


2.714 


2.736 


2.706 


r{Au6-Au7) 


2.839 


2.804 


2.832 


2.822 


r{Au7-Au8) 


2.772 


2.705 


2.729 


2.700 












r(Aui-Au2) 


2.828 


2.782 


2.822 




r(Au2-Au3) 


2.813 


2.782 


2.806 




r(Au3-Au4) 


2.789 


2.756 


2.787 




r(Au3-Aug) 


2.819 


2.774 


2.816 




r{Au5-Aue) 


2.886 


2.730 


2.759 




r(Au6-Au7) 


2.886 


2.852 


2.875 




r(Au7-Au8) 


2.724 


2.699 


2.728 














r{Aui-Au2) 


2.890 


2.829 


2.877 




r(Au2-Au3) 


2.739 


2.716 


2.740 




r(Au3-Au4) 


2.857 


2.811 


2.877 




r{Au3-Au6) 


2.927 


2.869 


2.910 




r(Au5-Au6) 


2.723 


2.706 


2.728 




r{Au6-Au7) 


3.008 


2.940 


2.979 




r{Au7-Au8) 


2.744 


2.717 


2.745 




Au^ 










r{Aui-Au2) 


2.850 


2.821 


2.849 




r{Au4-Au5) 


2.759 


2.740 


2.759 




r{Au5-Aus) 


2.853 


2.808 


2.849 




r(Au7-Au8) 


2.721 


2.707 


2.725 





Au^" 

r(Aui-Au2 
r(Aui-Au3 
r(Au2-Au7 
r(Au2-Aus 
r( Au4- Au5 
r{Au5-Au6 
Aug 

r(Aui-Au2 
r(Au2-Au4 
r(Au3-Au6 

r( Au3- Au^ 
r(Au4~Aufi 
r(Au4- Au7 
r(Au6- Au7 
r{Au6-Au8 

r( Aui - Au2 
r(Au2-Au3 
r{Au2-Au4 
r(Au3-Au4 
r(Au3-Au5 
r(Au3-Au6 
r{Au4-Au6 
r(Au4-Au7 
r(Au5-Au6 
r(Au6-Au7 
r(Au7-Au8 

r( Aui - Au2 

r(Au2-Au3 

r(Au2-Au6 

r(Au3-Au5 

r(Au4-Au5 

r(Au5-Au8 

r( Au7- Aug 

Au^^ + 

r(Aui-Au2' 

r(Au2-Au3 

r(Au2- Au4 

r(Au2- Au7 

r(Au3-Au8 

r(Au4-Au5 

r(Au4-Au7 

r(Au5-Au8 

r(Au7-Au8 



Brs(i/A 



2.817 


2.778 


2.812 


2.775 


2.658 


2.643 


2.664 


2.637 


2.978 


2.938 


2.968 


2.934 


2.817 


2.778 


2.812 


2.773 


2.658 


2.643 


2.665 


2.637 


2.784 


2.758 


2.784 


2.763 


2.774 


2.747 


2.776 




2.936 


2.887 


2.929 




2.920 


2.887 


2.913 




2.763 


2.727 


2.757 




3.()4{) 




3.t).-i8 




2.707 


2.690 


2.713 




2.796 


2.758 


2.791 




2.806 


2.779 


2.807 




2.732 


2.705 


2.736 




2.715 


2.697 


2.716 




2.840 


2.800 


2.844 




2.777 


2.745 


2.773 




2.779 


2.741 


2.773 




3.076 


2.943 


3.012 




3.378 


3.275 


3.287 




2.668 


2.654 


2.674 




2.733 


2.715 


2.740 




2.795 


2.753 


2.786 




2.962 


2.877 


2.933 




2.703 


2.681 


2.706 




2.886 


2.846 


2.875 




2.740 


2.723 


2.743 




2.970 


2.906 


2.960 




2.672 


2.670 


2.682 




3.033 


2.958 


3.011 




2.729 


2.711 


2.731 




2.696 


2.675 


2.702 




2.871 


2.844 


2.665 




2.936 


2.891 


2.936 




2.977 


2.934 


2.975 




2.786 


2.759 


2.786 




2.796 


2.765 


2.797 




2.931 


2.895 


2.930 




2.834 


2.799 


2.831 




2.767 


2.734 


2.759 
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Table 6: The bond lengths (in A) of the neutral clusters of nine gold atoms. 





B3LYP/A 


B3LYP/B 


B3LYP/C 


BP86/A 


r(Aui-Au2) 
r(Aui-Au3) 
r(Aui-Au4) 
r(Au3-Au4) 
r(Au4-Au5) 
r(Au4-Aug) 
r(Au6-Au7) 


2.719 
2.902 
3.161 
2.776 
2.769 
2.788 
2.868 


2.698 
2.850 
3.033 
2.745 
2.753 
2.746 
2.810 


2.726 
2.885 
3.092 
2.775 
2.774 
2.787 
2.851 


2.694 
2.857 
2.991 
2.747 
2.764 
2.748 
2.805 


Aui'CAi C4.) 
r(Aui-Au2) 
r(Au3-Au4) 
r(Au4-Au6) 


2.912 
2.922 
2.710 


2.856 
2.887 
2.688 


2.909 
2.919 
2.711 


2.851 
2.880 
2.684 


Aui^'CAiD^h) 

r(Aui-Au2) 

r(Aui-Au3) 

riAni-Aii/i) 

r(Au3-Au4) 
r(Au3-Au9) 
r(Au4-Au6) 
r(Au2-Au5) 
r(Au5-Au7) 
r(Au5-Au8) 


2.761 
2.956 
2.779 
2.752 
2.670 
2.891 
2.706 
2.950 
2.802 


2.733 
2.882 
2.751 
2.725 
2.658 
2.831 
2.687 
2.885 
2.759 


2.758 
2.930 
2.782 
2.756 
2.677 
2.868 
2.711 
2.930 
2.797 


2.737 
2.884 
2.753 
2.709 
2.650 
2.849 
2.676 
2.905 
2.755 


r(Aui-Au3) 
r(Au2-Au3) 
r(Aui-Au4) 
r(Aui-Aug) 
r(Aui-Au7) 
r(Au4-Au5) 
r(Au4-Au6) 
r(Au5-Au6) 


2.876 
2.745 
2.876 
2.844 
2.844 
2.729 
3.019 
2.679 


2.729 
2.729 
2.840 
2.755 
2.842 
2.753 
2.928 
2.645 


2.793 
2.746 
2.866 
2.834 
2.837 
2.726 
3.023 
2.686 


2.830 
2.703 
2.798 
2.796 
2.833 
2.674 
3.015 
2.675 



38 



This figure "Figlgold-JCP.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0412077vl 



This figure "Fig2gold-JCP.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0412077vl 



This figure "Fig3gold-JCP.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0412077vl 



This figure "Fig4gold-JCP.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0412077vl 



This figure "Fig5gold-JCP.gif" is available in "gif" format from: 



http://arXiv.org/ps/physics/0412077vl 



